SPECTRAL CHARACTERISATION OF AGEING: THE REM-LIKE 

TRAP MODEL 
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Abstract. We review the ageing phenomenon in the context of simplest trap model, 
Bouchaud's REM-like trap model from a spectral theoretic point of view. We show that 
the generator of the dynamics of this model can be diagonalised exactly. Using this 
result, we derive closed expressions for correlation functions in terms of complex con- 
tour integrals that permit an easy investigation into their large time asymptotics in the 
thermodynamic limit. We also give a 'grand canonical' representation of the model in 
terms of the Markov process on a Poisson point process . In this context we analyse the 
dynamics on various time scales. 

Key words: disordered systems, random dynamics, trap models, ageing, spectral prop- 
erties. 



1. Introduction 

The particular properties of the long term dynamics of many complex and/or disordered 
systems have been the subject of great interest in the physics, and, increasingly, the 
mathematics community. The key paradigm here is the notion of ageing, a notion that 
can be characterised in terms of scaling properties of suitable autocorrelation functions. 
Typically, ageing can be associated to the existence of infinitely many time-scales that 
are inherently relevant to the system. In that respect ageing systems are distinct from 
metastable systems which are characterised by a finite number of well separated time scales 
corresponding to the live times of different metastable states. 

Ageing systems are rather difficult to analyse, both numerically and analytically. Most 
analytical results, even on the heuristic level, concern either the Langevin dynamics of 
spherical mean field spin glasses, or trap models, a class of artificial Markov processes that 
in some way is to mimic the long term dynamics of highly disordered systems (see e.g. 

0)- 

One of the natural questions one is led to ask when being confronted with phenomena 
related to multiple time scales is whether and how they can be related to spectral properties. 
This relationship has been widely investigated in the context of Markov processes with 
metastable behaviour (see e.g. ^JJ El EES IZQ HU an d it would be rather interesting 
to obtain a spectral characterisation of ageing systems as well, at least in the context of 
Markov processes. To our knowledge, this problem has not been widely studied so far. 
The only paper dealing with the problem is the paper [23] by Butaud and Melin that have 
tackled one of the simplest trap models and on which we will comment below, and 18, and 
24 that investigate convergence to equilibrium in the Random Energy Model (REM). 

The present paper is intended to make a modest first step into this direction by analysing 
the relation between spectral properties and ageing rigorously in the REM-like trap model. 
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While this model may seem misleadingly simple, it has in the past provided valuable 
insights into the mechanisms of ageing, and it is our hope that the analysis presented here 
will provide useful guidelines for further investigations of more complicated models. 

The paper will be divided into two parts. In the first we analyse the REM-like trap 
model in the standard formulation of Bouchaud [Hj • In the second part we go a step further 
and reformulate the model in a slightly different way as a Markov process on a Poisson 
point process. This formulation makes the relation to the real REM more suggestive (see 
@] for a full analysis) , and allows in a natural way to study the dynamics of the model 
on different time scales. 

2. The REM-like trap model 

Let us begin to recall the definition of trap models as introduced by Bouchaud and 
DeanjH]. Let Q = (S,£) be a finite graph: S denotes the set of vertices and £ denotes the 
set of edges, E := {Ei, i G S} be a random field, called energy landscape and let Y(t) be 
a continuous-time random walk on Q with i?-depending transition rates q ,-: i.e. c$ »• > 
iff {i,j} G £ and 

F(Y(t + dt) = j | Y(t) = i) = a d dt. 

Setting n := Ylj^i c i,j an d Pi,j '■= (Hj/^u the random walk Y(t) can be described as 
follows: after reaching the site % the system waits an exponential time with expectation 
Tj and then it jumps to an adjacent site j with probability pij. In the trap model, the 
transition rates are assumed to satisfy the following properties: 

e Ei c id = e^c^ h V{i,j}e£, (2.1) 

Efa) = oo (2.2) 

where E denotes the expectation w.r.t. the random field E. Since in several physical 
experiments (see ^01) the system initially in equilibrium at high temperature T 3> T g is 
quickly cooled under the transition temperature T g and then its response to an external 
perturbation is measured, it is reasonable to consider Y(t) with uniform initial distribution. 
A classical time-time correlation function is given by 

U(t,t w ) :=P(Y(s) =Y(t w ), VsG [t w ,t w +t]) 

In order to observe ageing it is necessary to consider a thermodynamic limit, with the size 
of Q going to infinity, and possible a suitable time-rescaling. Rather recently, there have 
been a number of rigorous papers devoted to the analysis of trap models on the lattices Z 
[ElEniEl and Z d 0[12]. 

In this paper we consider the simples trap model, called the REM-like trap model (Hj 
that corresponds to choosing Q to be the complete graph on N vertices i.e. 

Gn = {Sn,£n), $n ■■= {1,2,. . . ,N}, £ N := {{i,j} i ± j G S N }, 

and to take as energy landscape a family E = {Ei : i G N} of independent exponential 
random variables with parameter a such that < a < 1. Given iV G N, let Yj^(t) be 
the continuous-time random walk on Qn with transition rates Cjj = e~ Ei /N for i ^ j. 
Setting Xi = e~ Ei the infinitesimal generator of the random walk is given by 



L 



N 



/ CAT-Dsi _xx _xi \ 

' N N N 

_X2 {N-1)X2 _X2 

N N ■ ■ ■ N 
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' N N ■ ■ ■ N 



(2.3) 
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The dynamics can be described as follows: after reaching the state i the system waits an 
exponential time of mean J\ e -E ' and then jumps with uniform probability to another 
state. Although strictly speaking the mean waiting time is given by j7zje Ei , we call 
Tj := x~ l = e Ei waiting time (the discrepancy is negligible in the thermodynamic limit 
TV Too). 

Note that Tj and xt have distributions respectively given by 

p(r)dr = OLT~ x ~ a dT (r > 1); p(x)dx = ax a ~ l dx (0 < x < 1), 
in particular E(rj) = oo. Moreover, the equilibrium measure is given by Heq(i) = 
r «/(SjLi T j) • We are interested in the out-of-equilibrium dynamic with uniform initial 
distribution. P^r denotes the law of this random walk given a realization of the random 
variables E{. 

Ageing in the REM-like trap model is manifest from the asymptotic behaviour of the 
time-time correlation function 

U N (t,t w ) :=¥ N (Y N (s) =Y N (t w ), VsE [t w ,t w + t]) (2.4) 

Namely, as shown in j^j, for almost all E and for all 9 > 

lim lim U N (9t w ,t w ) = f u^il-u^du (2.5) 

t w '\ooN]oo IT J 8 

1+0 

Our main aim here is to show that the ageing behaviour of the system, derived in 
[2] using renewal arguments, can be obtained solely from spectral information about the 
generator Ljv- The method developed below will allow us to get further information on 
Yff(t) from the spectral properties of L^r. In particular, given a function h on (0, oo), it 
is possible to describe the asymptotic behaviour of EN(h(x]y(t)) and E7v(/i(rAr(t)), where 
Ejv denotes the expectation w.r.t. Pat and xn (t) , tn (t) are defined as 

x N (t) = x k , r N (t) = r k if Y N (t) = k. 

These results will allow us to investigate the property that the system with high probability 
visits deeper and deeper traps, i.e. sites with larger and larger waiting time Tj (see 
subsection 12. 2 j) . 

We start by giving a complete description of the eigenvalues and eigenvectors of Ljv. 
Let fi = fiN be the measure on Sn with fi(i) = x^ 1 = Ti. Note that ~Ln is a symmetric 
operator on L 2 (fi) and trivially LtvI = where I is the vector with all entries equal to 1. 
The following proposition is based on elementary linear algebra: 

Proposition 2.1. Let x%, X2, ■ ■ ■ ,xn be all distinct. Then, L^v has N positive simple 
eigenvalues = Ai < A2 • • • < Xn such that 

{A 1 ,A 2 ,...,A iV } = {AGC : 0(A) = 0} 

where 4>(\) is the meromorphic function 

N A 
3=1 Xj A 

// the Xi are labelled such that x\ < X2 < ■ ■ ■ < xn, then Xi < Aj+i < Xj+i for i = 2, . . . , N . 
Moreover, for any i = 1, . . . , N , the vector ipw G M. N defined as 

^ ): =—^T> forj = l,...,N 



1 



ANTON BOVIER AND ALESSANDRA FAGGIONATO 



is an eigenvector of L^v with eigenvalue Aj. ip^\ ■ ■ ■ , ifi( N ' form an orthogonal basis of 

Since the Xi have a absolutely continuous distribution, we trivially have the 

Corollary 2.2. The assertions of Proposition HOI hold with probability one for all N. 

Proof. Let A be a generic eigenvalue and let us write a related eigenvector vj as ip = 
a(l, . . . , 1)* + w where Y2j=i w j = 0- Since (LjvV')j = x j w ji we have to solve the system 

XjWj = Xa + Xvjj, Vj = 1, . . . , N. (2.7) 

Since x\, . . . , xn are distinct, it must be true that a/0 (otherwise we get ip = 0). Without 
loss of generality, we set a = 1. Note that A ^ Xj for j = 1 . . . TV, as otherwise (|2.7|) would 



imply A = = Xj. Therefore we get Wj = x x _ x - Since it must be true that Ylf=i w j = 0> 
we get that A is an eigenvalue and ifi, with yjj = x \ , the corresponding eigenvector, 
iff 0(A) = 0. This imply that <f> has at most N zeros. Since 0(0) = 0, and, for real A, 
limAj.Xi 0(A) = — oo, limAfxi 0(A) = °°> we get that has exactly N zeros. The conclusion 
of the proof is trivial. □ 

Proposition 12.11 has the following simple corollary. 

Corollary 2.3. With probability one, the spectral distribution o~n := Kv I ^ =l 5\ j converges 
weakly to the measure ax a ~ 1 dx on [0, 1]. 

Remark. The results of Proposition 12.11 are incompatible with the predictions of |25j . 
The discrepancy is particularly pronounced in the case of the eigenf unction. The reason 
for this is an inappropriate use of perturbation expansion in |25| . We will explain this in 
some detail in an appendix. 

We will now show that Proposition 12.11 allows to derive the asymptotics of the auto- 
correlation functions easily. In fact, it contains far more information on the long time 
behaviour of the systems some of which we will bring to light later. 

Recall that pt(i,j), the probability to jump from i to j in an interval of time t, can be ex- 
pressed as Pt(i,j) = (e _tLjv ) . .. In particular, by writing vt for the probability distribution 

of Yjsr(t) and thinking of the Radon derivative ^7 as column vector, 



*1 = c -fLy dl/ 

d/j, dfj, ' 

thus implying 

N ^ duo ,/,(fc) ^ 

P- = y t: e-*v° (2.8) 

dfi ^<^ fc \V (fc) > 

The above formulas are true for an arbitrary initial distribution. Taking uq the uniform 
distribution, by Proposition 12. II we get 

, N N 

dVn 



'O 



Y,lk4 k \ where 7fc " 1 :=< ^ k \^ >= ^ 
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Then, by Proposition EU and l|2~%]) . 



^ —Xl-ty 



^»(M.) = EE 1: ^ : «- ¥ " , (") 

j=l k=l J 

N N -Afct 

E N (h(x N (t))) = ^2^1^-^h( Xj ) (2.10) 

j=l k=l J 

The above formulas (that may appear rather ugly at first sight) admit a nice complex 
integral representation through the following lemma. 

Lemma 2.4. Let 7 be a positive oriented loop on C containing in its interior Ai, . . . , Ajv- 
Let g be an holomorphic function on a domain D C C with 7 C D. Then, for any 
j = l...,N, 

^x^-Afc 27rz 7 7 <j)(\)(xj - A) v 7 

Proof. Let us set X := {xi, . . . , xn} and A := {Ai, A2, • • • , Ajv}- Then 0(A) is an holo- 
morphic function on C \ X, where <f>'(\) = YljLi ( x -a)S > an< ^ ™ particular = 7J" 1 . 
Moreover, the function [<p{\) (xj — A)]" 1 that is a priori defined on C \ (X U A) can be an- 
alytically continued to X to a meromorphic function with simple poles only at the points 
of A. Now the conclusion follows from a trivial application of the residue theorem. □ 

We can obviously use Lemma 12.41 to rewrite formulas (|2.9|) and ()2.1(jj) in the form 

where Av,- denotes the average over j = 1, 2, . . . , N . 

The above integral representations of H^{t,t w ) and ~K^{h{xt)) have two main advan- 
tages. First, the appearance of averages allows to compute their limiting behaviour as 
N j 00 easily by using the ergodicity of the random field E_. Second, by means of the 
residue theorem, their Laplace transform can be easily computed in order to derive the 
asymptotic behaviour of U]\f(t,t w ) and E7v(/i(a^)) for N,t w ,t 3> 1 (see subsections 12.11 

E2>. 

A much more general derivation of the above integral representations is discussed in 
Appendix O 

2.1. Ageing behaviour of ILv(i, t w ). 
Proposition 2.5. Let us define 

U(t,t w ) : = — / — ) x - x ( d\. (2.14) 

where K x is the expectation w.r.t. the measure ax a ~ 1 dx on [0,1] and 7 is any positive 
oriented complex loop around the interval [0, 1] . Then 

lim IIjv(t, t w ) = II(t, t w ) yt,t w , a.s.. (2-15) 

Art 00 
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Proof. Recall (|2.12j) and fix < 5 < 1/2. Due to analyticity, we can choose the integration 
contour 7 to have distance 1 from the segment [0, 1]. For each A G 7, the random variables 
(xj — A) -1 , j G N, are i.d.d. and bounded. Therefore, for a suitable positive constant 
O0, 



Av 



N 

j=1 1" • - A 



E, 



(- ) > N~2+ d ) < e 

A — rf -' 



VA G 7. 



(2.16) 



Since for each x G [0, 1] and A G 7, — A) _1 | < 1, a simple chaining argument allows to 
deduce from the pointwise estimate (|2.1fij) uniform control in A. With the Borel-Cantelli 
lemma one can then infer that, a.s., 



sup 

A£7 



Av 



N 



1 



1 



E *(— ) 



"J 



VA G N. 



(2.17) 



Similar arguments show that, a.s., given M G N there exists a constant cm such that 

N Z N 



sup sup 

M-l<t<M AG7 



Av 



- A 



e 4t— -) 



-X-Xj' 



<c M N 2 



VA^ G N. 



2.18) 



Note that, for each A G 7, Av^ 1 (x J - — A) -1 is a convex combination of points of modulus 
larger or equal than 1/2, contained in a angular sector with angle non larger than a suitable 
constant c < ir. In particular, |Av^ 1 (x J — A) _1 | > d > for all N. From here the assertion 
of the proposition follows from Lebesgue's Dominated Convergence Theorem. 



□ 



Given 9 > we are interested in the limit of H(9t w ,t w ) as t w 1 00. This will be done 
using the Laplace transform of IL(9t w ,t w ), 



11(6, u) :-- 



-LOt v 



n(9t w ,t w )dt w , (»(w)>0). 



The computation of this Laplace transform is trivial if we use the integral expression 
(EH- 




Let ijj G C with 5R(cj) > and fix a positive oriented loop 7 around the segment [0, 1], 
such that 7 C {z G C : 9E( z) > -SR(cj)} as in fig. O Then, R(w + A + x9) > 0, for 
x G [0,1] and A G 7, so that (|2.14j) implies 



2ni 



11(9, io) = E x (— A(A - x)(X + w + flar)Es( 



A — x' 



dA 



Here Ej, and E s denote the expectation w.r.t. to the measure ax a ~ 1 dx on [0, 1]. 

Let us consider the change of variable z = \ and write 7 for the path 7 with inverted 
orientation (i.e., positive oriented w.r.t. A = 00). Then we get 



11(9, oj) =^(2^ / ( I- - '■'')( 1 + + l9.r)Ej 



1 — zx' 



dz). 
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Given x € [0, 1], the integrand is a meromorphic function in C \ [1, oo) that has only a 
single pole of order 1 inside 7, namely at z = — (u + x0)~ l . By the residue theorem we 
get 

Lemma 2.6. The r.h.s. of h2.iy\) is well defined and holomorphic for any uj G C\(— 00, 0]. 
In particular, the function 11(0, uj), defined for dt(u) > 0, can be analytically continued to 
the set C \ (—00, 0]. 

Proof. As proved in [TS|, Chapter 3, the Laplace transform 11(0, uj) is holomorphic on the 
set of convergence points. Therefore, we only need to show that the r.h.s. of ()2.19j) is 
well defined and holomorphic on $S{u) 7^ 0. Let us assume ^s(ui) > a > 0. Then, trivially, 
Vx,x G [0,1], 



CO + XV + X 



eB:={zeC : z= \z\e lti with < 6 < 6 , \z\ > c} 



for suitable constants c, #0 depending on a and such that #0 < ir. Moreover, since 
< cUa) < I u + * d _ | < C2 ( a ) Va > 0,Va; : 3f(cj) > a (2.20) 

' LU + X0 + X 

By (|2.20|) and the geometry of B, we have that K x ( rjr^frg ) is well defined and has distance 
03(a) > from the origin. Moreover, |cj + x8 + x| > 3?(o>). Therefore, the r.h.s. of (|2.19)1 
is well defined and, due to the previous estimates and Lebesgue's dominated convergence 
theorem, it is continuous on {3?(w) 7^ 0}, thus implying continuity on C \ (—00, 0]. 

We recall Morera's theorem: if f(oj) is defined and continuous in a open set O C C 
and if J f dw = for all closed curves 7 in Vt, then f(u>) is holomorphic in S7. Therefore, 

using Fubini's and Morera's theorems, one can prove that the function ~E X ( J^q^ ) is 
holomorphic on C \ (— 00, 0]. The proof can be concluded by a second application of the 
same theorems. □ 

In what follows, we keep the notation 11(8, to) for the analytic continuation of the Laplace 
transform. The next lemma describe the behaviour of H(9, to) near the origin. Using the 
Laplace inversion formula, we then derive from this result the asymptotic behaviour of 
11(0,4,) as t u | 00. 

Lemma 2.7. For any 6 > let us set 



Moreover, let us define 



vr J_e_ 

e+i 



A:={re^ : r > 0, \<f>\ < ^vr}. (2.21) 



Then, for a suitable positive constant c > 

|n(6>,a>)| < c\uj\~ 1 , Vw G A : \u>\ > 1 (2.22) 

\il(e,u) -A(9)/u\ < c\uj\- a , VweA:\u\<l. (2.23) 
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Proof. The first estimate (j2~2~2jl follows trivially from (|2~T3j) and Let us prove (|2~2T?|) 

for u; E .A and |cj| < 1. 

In what follows, cq,ci, . . . denote some suitable positive constants depending only on 
0. Moreover, given zGCwe denote by J Q Z and Jj 30 the integrals over the paths {sz : < 
s < 1} and {sz : s > 1}, respectively. We extend the functions z~ a and defined on 

(0,oo), to C \ (—oo,0] by analytic continuation. Then (|2.19|) implies 



Let us define 



1/uj 



X 



a-1 



1/uj 



X 



a-1 



[l + x(l + 0)] [l+x9] 



[l + x(l + 0)] [l+x9] c 



dy l dx 

1 + xO + y y J 

"(i+*e) \ -1 

ay ax. 



(2.24) 



o 



1 + y 



B := i) EC s.t. u) £ A, x = — for some s : < s < 1}. 

w 

Since (u>(l + £0)) 1 G A D {2 : |z| > Co}, we obtain 

„ct-l 



1 + 1/ 

a-1 



dy 



> ci, 



_^ 1+1/ 



< C2 |u>(l + x£ 



ll-Q 



Let B{9) be defined as 
5(0) := 



OO yOl-1 



1 + y 



dy 



TT 



(2.25) 
(2.26) 

(2.27) 



u- a (l - u) a ^du = 

sin(7ra) 

(note that the above second identity follows from the change of variable u = y(l + y) -1 , 
while the last one is well known in the theory of the Gamma function) . By means of ()2.25|) 
and (|2.26|) we obtain 



1 



x a 1 dx 



< 



B(9)J (1 + 1(1 + 0)) (l + x9) a 

l^l"" 1 d\x\ 



1-a 



(2.28) 



Since 



\l + x(l + 6)\\l + x0 

x a ~ 1 dx 



2a- 



1-a 



< C4 \u)\ 



f i (l+x(l + 0)) {l + x9) a 
and, using analyticity and integrability of the singularities around z = and z = 00, 

x a ~ l dx _ f°° x a ~ l dx 

(1 + x(l + 9)) (1 + x9) a J (l + x(l + 9)) (l + x9) a, 

we get 



/ 

J 50. 



ujU(9,tt) 



1 



„a-l 



-dx 



< c 5 |w 



1-a 



B{9) J (l + x(l + 0))(l + x0) a 
Using the chance of variables v = x^ 1 + 9 and u = v(l + v) , we obtain 



.1" 



a-1 



(l + x(l + 0))(l + x6 



-dx 



u -«(l _ u^du 



e+i 
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which implies the assertion of the lemma. □ 

Lemma 12.71 and Proposition I A. II allow us to conclude the proof of the ageing behaviour 
of Il N (t,t w ): 

Proposition 2.8. For almost all energy landscapes E_, given 9 > 

lim lim U N (9t w ,t w ) = sm ( 7rQ ) f 1 u - a {l-u) a ~ l du. (2.29) 

t w ]ooN]oo TT I e 

i+0 

2.2. Visiting deeper and deeper traps. 

In this section we use the integral representation (|2.13|) in order to study the proba- 
bility that the system at time t is in a deep trap, i.e. in a state with large waiting time. 
In Proposition 12.91 we first prove that the probability to be in a site with waiting time 
smaller than O(l) decays as thus implying the ageing behaviour of other correlation 
functions described in subsection |3 In the second part, we will investigate the random 
variable tx]y(t) and show that, for almost all E, it has a weak limit as N f oo and then 
t | oo. As consequence, with high probability at time t the system is in a state of waiting 
time 0(t) as stated in Proposition 12. 1UI 



Reasoning as in the proof of Proposition (|2.5|) . we can prove, for almost all energy 
landscapes E, that, given a function h on [0, 1] that can be uniformly approximated by 
piecewise C 1 functions, 

H(t) := lim E N (h(x N (t)) = — / - ^ ~y~ — dX Vt > 0, (2.30) 

where 7 is a positive oriented loop around [0, 1]. 

Since H(t) is a bounded function, the Laplace integral H(uS) := / °° H(t)e~ wt dt is 
absolutely convergent when > 0. By the same arguments we used to derive ()2.19j) . 

it is simple to deduce from the integral representation (|2.30|) that 

H(u) = - J -\^ (2.31) 

UJ f 1 T a-lrj r 
JO uj+x x ux 

In the following Proposition we concentrate on the case h{x) := I x >s- By (|2.31j) . we can 
give precise information on the asymptotic behaviour of the probability to be at time t in 
a site with waiting time smaller than 1/5: 



Proposition 2.9. Let 



■= it =5=i ; , c(a) := I y^e^dy (2.32) 



roo 

JO 1+x ux 

Then, for almost all energy landscapes E, 

lim s 1 "" lim F N (x N (s) > 5) = B{5)/c{a). (2.33) 

sjoo N^oo 

Finally, we show that with high probability at time t the system is in a trap of depth of 
order 0{i). In particular, the random variables tx^(t) converge weakly to a nonnegative 
random variable as N f 00 and then 1 1 00 a.s.. This result corresponds to the convergence 
of expectation of bounded continuous functions and due to Lemma 12 . 1 1 1 such a convergence 
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can be extended to the larger class of bounded piecewise continuous functions, which is 
more suitable in order to investigate the phenomenon of visiting deeper and deeper traps: 

Proposition 2.10. Let Z be the only random variable with range in (0, oo) having Laplace 
transform 

ez sinM /- 1 u - a{1 _ ur -x du 

Then, for almost all energy landscape E, given a bounded piecewise continuous function h 
on (0, oo), 



tx rl Hxt) x a-l dx 

lim lim E N (h(tx N (t))) = lim — / J ° \ x — dX = E(h(Z)) (2.34) 



In particular, for almost all energy landscapes E, 

lim lim P(^^ > u) = F(Z < u" 1 ), Vu > 0. 

ttoo 7V|oo t 

Proof of Provosition \2.fA We have to prove that lim^oo s 1 ~ a H(s) = B(5)/c(a) where H 
is given by (|2.30f) with h{x) := I x >5- As in the proof of Lemma 12.61 we can show that the 
r.h.s. of (|2.31|) is well defined and holomorphic on C\ (— oo,0]. We keep the notation H 
for this extended function. By the change of variables x = toy, we get 

1 „q-1 i- „.a— l 

dx = LO a ~ 1 / f dy (2.35) 



where 7^ is the oriented path {s/u;}o<s<i- Let % be the path {s/u>} s >q. By analyticity 
and integrability of the singularities at z = 0, z = 00, we have 

yOt— 1 POO yOi—\ 

= / dy 

/%, 1 + 2/ V 1 + y 

Let us define .A := {re lS : < r < 00, |#| < jTt}. Then, for a suitable constant c\, 

yOL-l 

dy\ < ci\uj\ l ~ a , VujgA: \u\ < 1, 



7w 



1 rpQ — 1 / f 00 1 



implying 

f 1 ^ dx = UJ ^l I ^ dy + o(\u\^ a )), (2.36) 

where A = i? + 0(1/N) is understood to mean that there exists C < 00 such that 
\A - B\ < C/N. Trivially, 

,1 x^ dx = , 1+ , r 1 xQ - 2cfa (2>37) 

Let us note that the estimate of error terms in (|2.36|) and in (|2.37|) is uniform in u E A, 
\u\ < 1. Then, from (TOTT) . (t2~3Hl) . (t2~371) we get 

|a; a #(u;) - B(5)/c(a)\ < c 2 \uj\ 1 - a , VuGi: |w| < 1. (2.38) 

Since trivially |.ff(u;)| < C3 1 u-' | — 1 for ui E A with \u\ > 1, the assertion of the proposition 
follows from Proposition IA. 11 □ 
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Proof of Proposition \2.l(K As discussed before (|2.H()|) . one can show, for almost all 
energy landscapes E, that, given a piecewise continuous function h on (0,oo), 

i r -tx f 1 h^l x a-i dx 
$ t (h) := lim E N (h(tx N (t))) = — / - Jo rl x - x — dX Vt > 0, 



where 7 is a positive oriented loop around [0, 1] . Note that $f defines a positive linear 
functional on the space of continuous functions on (0, 00), decaying at 00 and satisfying 
$t(l) = 1. Therefore, the Riesz-Markov representation theorem (see Theorem IV. 18 in 
26 ) implies that $t(/t) = m(h), for a unique Borel probability measure fit on [0, 00). In 
particular, there exists a random variable Z-t on (0, 00) such that 

lim txpfit) — ► Zt weakly , Vt > a.s.. 

If we take h(t) = e~ w , then $ t (h) = m(h) = U(9t,t), with II defined as in (|2~Ti|) . That 
means that H(6t,t) is the Laplace transform of Zp As proved in subsection 12 .!( 

limn(gf,t) = sm ^ f 1 v r a (l-u) a - 1 du:= f{9). 

We state that f(6) is the Laplace transform of a random variable Z with range in [0, 00). 
To this aim we apply the criterion given by Theorem 1, Section XIII. 4 in By (|2.27|) . 
/(0) = 1. Moreover, f <yl \9) = — sm ^ a ^ Q-°(i -|- thus implying (by trivial computa- 
tions) that (—l) k f( k \6) > 0. This completes the proof of our statement. 

Since the Laplace transform of Zt converges to the Laplace transform of Z as t goes to 
00, we have that Zt converges weakly to Z, thus implying the limit ()2.34[) whenever h is 
a bounded continuous functions on (0, 00). Finally, due to Theorem 5.2 in [7j, the limit 
remains valid if h is a bounded measurable function whose set of discontinuity points has 
zero measure w.r.t. the distribution of Z. Therefore, Lemma 12.111 allows to prove ()2.34|) 
for h bounded and piecewise continuous. □ 

Lemma 2.11. The distribution function F(z) := ¥(Z < z) of the positive random variable 
Z is continuous. 

Proof. Trivially F is increasing and right continuous. Therefore, it has a countable set of 
points of discontinuity. Moreover, by the Laplace inversion formula (see ^Jj, XIII. 4), if x 
is a point of continuity, then 

(-ay 



F{x)= lim £ ^/W( fl ). 



n<aa; 



Given s = 0,1,2,... and 7 > let 0,(7) > be such that D s a a 7 = (— l) s c s (j)a 7 ~ 
Then the Leibniz formula implies 

n 

(-l) n ^(a- Q (l + a)- 1 ) = Y,c s (a)c n ^ s (l)a- a - s (l + a)- 1 - n+s < {-l) n Dy-^ 1 . 



s=0 



Since /W(a) = — sm ^ ra ) Q Q (i -|- a ) 1 ) the above estimate implies 



n— 1 



(-l) n / (n) (a) = |/ (n) (a)l < + a ) e ~ n ~ a ' Vn > 2 - 



k=l 
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In particular, given two points of continuity < x < z we have 

n— 1 



F(z) -F(x) < limsupa-" ^ IjJ(i + ^). ( 2 - 39 ) 



a— \ n k 

ax<n<az k=l 



One can prove that the sequence 11^=1 e n (-*- + §) ^ s convergent (see [Q, Chapter 5, Section 
2.4). Let us denote its limit by c a and let 7 be Euler's constant 

/ll 1 , \ 
7 := lim 1 + - + - H 1 log n . 

ntoo V 2 3 n / 



Then we can write 

- ff(l + r) = e-^+H-'+^M-D) (J^T TT e _ f (1 + a 

n - L - L fc n - LJ - A; 

fc=i fc=i 

In particular, in (|2.4()j) we can substitute Jlfc=i e "(l + f ) with c Q with an error term in 
(t2~331) bounded by 

n— 1 n— 1 

const. a~ a (az — ax){ax)~ l+a | e~^(l + ^) — c a < c(x, z) e~^(l + ^) — c 

fc=l ' fc=i 

which is negligible as a | 00. Therefore 

F(z) - F(x) < lim sup c Q a- a e 7Q V (n - < c'(z Q - x a ) (2.41) 

a^oo , 

ax<n<az 

for a suitable positive constant d . Since (|2.41() is valid almost everywhere and F is 
monotonic, we have that F is continuous. □ 

3. Other correlation functions 

In this section we study the asymptotic behaviour of some other time-time correlation 
functions U^\t,t w ), H$(t,t w ) for which deep traps play a special role. This section is 
mainly a preparation of what is to follow in the second part of the paper. 

Given 5 > 0, we define the set of sites with small waiting time as Dn '■= {i '■ Xi > 
5, i = 1, . . . , N}. Moreover, we set 

U^\t,t w ):=F N (Y N (u)eD N \/u e {t w ,t w + t] s.t. Y N (u) / Y N (u~) ) (3.1) 

n% ] (t,t w ) :=F N (Y N (u) e D N U{Y N (t w )} Vu € (t w , t w + 1] s.t. ^ x N (vT) ). 

(3.2) 

Given a subset A C 5tv, i G i and s > 0, let ifN,A(h s ) be defined as 
<PNAh*) :=Vn{Y n (u) eAVue [0,s] \Y N (0) = i), 

then 

N ,.f —sx - 

n§\t,t w ) = IL N (t,t w ) + J2 F N(Y N (t w )=j) ds Xj& N ' <PN,D N {i,t-s), (3.3) 

j=i Jo ieD N 

N 

n$ (M™) = ^2¥ N (x N {t w ) = xj) ipN,D N u{j}{j,t) (3-4) 

where the first identity can be derived by conditioning on the first jump performed after the 
waiting time t w and by recalling the following realization of the dynamics: after arriving 



SPECTRAL CHARACTERISATION OF AGEING 



13 



at the state i, the system waits an exponential time with parameter Xj and after that it 
jumps to a site in 5/v with uniform probability. 

The following proposition is mainly a consequence of the phenomenon of visiting deep 
traps with higher and higher probability. To this aim recall Proposition 12.91 or simply, as 
a consequence of Proposition 12.101 that 

lim lim F N (x N {t) > e) = 0, Ve > 0. (3.5) 

tf oo JVfoo 



Proposition 3.1. For almost all x 



lim sup|II$(M m ) - U N (t,t w )\ = fori = 1,2 (3.6) 



twloo £>() 

Proof. We consider first the case i = 1. 
We claim that for any u > and i E -Dtv, 



<PN,D N (i,u) < exp(-5u(l - -^r-))- ( 3 - 7 ) 

In order to prove such a bound, we introduce a new random walk Y^(t) having generator 
L* defined as the r.h.s. of (|2.3|) . where X{ is replaced by 5 if i G D^. By a simple coupling 
argument one gets 

<PN,D N (i,u) < <p* N)DN {i,u) 
where the function f* ND is the analogous of ipN,D N (i,u) for the random walk Yjy(t). At 
this point, it is enough to observe that (p* N Dn equals the r.h.s. of (|3.7j) . 
Now fix e > 0. Then, due to (jH3J) and $77) . 



\U { ^\t,t w )-U N (t,U 



<¥ N {x N {t w )>e)+ W N (Y N (t w ) = j)^ J2 [ fN,D N (i,t-s)ds 



j : Xj <e i£D N 

< ¥ N {x N {t w ) > e) + e I e" 5 "( 1 ~^ d )dn. 
Jo 

By the law of large numbers, 



(3.8) 



lim / e v n ) < cxd, a.s.. 
AT T oo J 

The proposition now follows from the fact that e is arbitrary and from (|3.5j) . 

To deal with case (ii), one proceeds in essentially the same way, decomposing the path 
of the process at its returns to the point Xi, and summing over the number of these returns. 
One finds easily that the case when the process does not leave X{ for the entire period t 
dominates, leading to the assertion of the proposition. We leave the details to the reader. 

□ 

4. The REM-like trap model on a Poisson point process. 

In this section we consider a slightly different formulation of the REM like trap model 
that betrays more directly its connection to the REM dynamics (see E]) and that 
offers a somewhat more natural insight in the role of time scales in the analysis of ageing 
systems. Let us consider a Poisson point process V = ^ • 5g 4 on M with intensity measure 
ae~ aE dE, where < a < 1. Note that such processes arise naturally as the extremal 
process of sequences of random variables. Trivially, a.s. the support of V is an infinite 
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set of points, whose maximum is almost surely bounded from above. Thus we can label 
the points in the support of V in decreasing order: E\ > E% > . . . . Then, the energy 
landscape E is defined as E = (E\,E2, ■ ■ ■)■ We want to define a random process on the 
support of this point process that jumps "uniformly" from any point to any other point in 
the support. To do this, we need to introduce a cut-off. Here we fix an energy threshold 
E and set 

N E = max{z : E { > E}. 

Note that Ne is a Poisson random variable with expectation e~ aE . Moreover, the proba- 
bility that Ne > 1 can be made as small as desired when E is chosen small enough, as we 
assume in what follows. 

Let Qe = (<Se,£e) be the graph with 

S E := {1,2,... ,N E }, £ E :={{i,j} : i^jeS E }. 

Since here we want to investigate the effect of time rescaling, we introduce a time unit 
ro = e E °. Then, the continuous-time random walk Ye(£) is the random walk on Qe 
having uniform initial distribution and such that, after arriving at site i G £e, it waits an 
exponential time with mean j^^je Ei / tq and then jumps with uniform probability to a 
different site of Se- In particular, the Markov generator for the above defined random 
walk is given by in (|2.M|) with N := Ne and Xi := T§e~ Ei (since for E <C 0, ^e-x ~ 1 
when referring to waiting time we disregard the coefficient j^ti as i n section |2j). Although 
Ye(£) depends on to, our notation does not refer to such a dependence. In what follows 
we denote by the probability measure on the paths space determined by Ye(-), and by 
E# the corresponding expectation. 

Note that the physical waiting time (the absolute one) for the system at state i is given 
by Tj := e Ei while in the above dynamics the waiting time is Tj := Tj/ro, thus in agreement 
with the choice to consider tq as our new time unit. In what follows we consider, when 
taking the thermodynamic limit E j — oo, three different kinds of time rescaling: To fixed, 
r := e E (that is E = E) and t j after E [ — oo. 

As in section 12 we are interested in the asymptotic behaviour of time-time correlation 
functions. In particular, let us introduce here the correlation function 

U E {t,t w ) :=Fe{Ye(s) =Y E (t w ), VsG [t w ,t w + t]). 

We will prove that when To is fixed the system exhibits fast relaxation, thus excluding 
ageing behaviour (see Proposition 14. 2|) . At the other extreme, the scaling To = e E cor- 
responds to the implicit choice made in the standard Bouchaud model considered in the 
previous sections. In fact with this choice the system can be thought of as a grand canoni- 
cal version of the original REM-like trap model and all the results of the previous sections 
carry over. Finally, we consider the third scaling: to j after E J, — oo. In Proposition 14.61 
we show that when performing such limits the correlation function Il£:(t, t w ) converges to 
f(6) where 9 = t/t w and f(9) denotes the r.h.s. of identity (|2.'29jl . that is the limiting 
behaviour of the correlation function TlE(t,t w ) is trivial. At this point a simple consider- 
ation is fundamental. If we assume that the physical instruments in the laboratory have 
sensibility up to the time unit To, then it is natural to disregard jumps into states with 
physical waiting time Tj = e Ei much smaller than To. Therefore, a time-time correlation 
function much more interesting than ~n.E(t,t w ) is the following one, where 5 > is fixed 

U { ^(t,t w )=F E {xE(u)>5 Vue(t w ,t w + t]: x E {u) ^ x E (u~) ). 
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where x E {t) := x k whenever Y E (t) = k. In Section [5] we prove that Tl^\t,t w ) exhibits 

ageing behaviour: 11^/ (Ot w ,t w ) converges to the above f(9) after taking the (ordered) 
limits E I — oo, tq j and t w ] oo. 

Finally, in this section we discuss the asymptotic spectral behaviour for the above time 
rescalings. We will show that ageing appears whenever the limiting spectral density has a 
singularity of order 0(x a ~ 1 ) at 0. 

Let us recall some properties of the Ppp ^ 5 Xi with intensity measure aTQ°'x a ~ 1 dx 
on (0, oo) which will be frequently used below. Given M > the truncated Ppp 
Y2x <M $3>i can k e realized as follows. Let hm be a Poisson variable with expectation 
= Jq 1 aTQ°x a ~ 1 dx and let Xj, i G N, be i.i.d. random variables on [0, M] with 
probability distribution p(X)dX = aM^X^dX. Then 

(4-1) 

Xi <M i=l 

in the sense that the above point processes have the same distribution. In particular, by 
taking M = T§e~ E , we get 

n* 

£ 4* (4-2) 

i<N E i=l 

where n* E is a Poisson variable with expectation e~ aE and G N, are i.i.d. random 
variables on [0, T$e~ E ] with probability distribution p(X)dX = e aE aTQ C1 X a ~ 1 dX. 
Notation It is convenient to introduce the random walks x E (t), T E (t) defined as 

x E {t) := x k , T E (t) := T k if Y E {t) = k. 

We denote by j E the positive oriented loop having support 

supp(7£) = {x±i : x G [-1, T e~ E }} U {-1 + bi : \b\ < 1} U {r Q e~ E + 1 + bi : \b\ < 1}. 

Moreover, we call the infinite open path, oriented from 00 + i to 00 — i, having support 

supp(7oo) = {x ± i : x> -1} U {-1 + bi : \b\ < 1} 

( E) ( E) ( E) 

Finally, for given E, = X\ < X 2 < ■ ■ ■ < X N ^ , are the N E distinct eigenvalues of 
the infinitesimal generator ~L E (see Proposition [2^). 

4.1. To fixed. 

Let us first observe that ^4=1 T i < 00 for almost all E_. In fact, since the Ppp Yli &n 
has intensity measure aTQT~( 1+a >dT on (0,oo), 



E(|{i: n>l}|) 



/"OO f'\ 

/ aT Q T- (1+Q W < 00, E( V n) = / ar^T^dr < 00. 
Jl i:n<l J ° 



Whenever ^£1 Tj < 00, it is simple to derive the asymptotic spectral behaviour of the 
system from Proposition 12. II and to show its fast relaxation, thus implying the absence of 
ageing: 

Proposition 4.1. For almost all E_, 

Ne 00 

lim y 5 (e) = / 5\ . vaguely in Ai([0, 00)), (4.3) 

E ^U x u 
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where «M([0, oo)) denotes the space of locally bounded measure on [0,oo) and 

A 

A - x k 



oo , 

{0 = Ai < A 2 < A 3 < . . . } = {A e C : y = 0} 



k=l 

Proof. In what follows we assume that T « < °°i which is true a.s.. Then the function 
=i Xk'—X ls wei l defined on C \ {xj : % > 1} and has non negative zeros 
= Ai < A2 < . . . , such that Xj_i < Aj < Xj for any i > 1. At this point it is enough to 
show that 

lim Af } = Aj, Vi = l,2,.... 

£|-oo 

The assertion is trivial for £ = 1. Suppose that iVg > i > 1 and set ^(A) := Sfe=i x fc -A ' 

Due to Proposition 12.11 A| is the unique zero of i/jeW i n the interval (xj_i,Xj). In 
particular, 

^B(Ai) = M^i) ~ ^(Af ) ) = P ^(A)dA. 

JA^ 

Since V>e(A) > ^tz^— ^3 for all A G (a^-i,^), we get 

|Af 5 -Ai| < (si-Xi_i) 2 |^(Ai)| 

and therefore the assertion follows by observing that the identity </>oo(\) = implies 

00 ^ 

|<MAi)|< V |0 as £ I -00. 



X L> X i) 

k=N E +l 



□ 



lim lim F E (x E (t) = X j) = =^ Vj = 1,2,... (4.4) 



Proposition 4.2. For almost all E_, 

lim lim P E (x E (t, .. , , 

thus implying 

lim lim n E (0i W) i w ) = V6> > 0, (4.5) 

t-w^oo El— oo 

lim lim n g (t,t w ) = Vt>0. (4.6) 

tmfoo El— 00 /_^i=l T i 

Proof. In what follows we assume that E satisfies Yli T i < 00 • By setting /i(x) = H a .= a ;- in 
(|2.13|) we get the integral representation 

PeM) = Xj) = (4.7) 

By applying the residue theorem (see the arguments used in order to derive (|2.19)) ). it is 
simple to compute the Laplace transform Fe{oj) = Jq°° P.E(x.E(i) = Xj)e~ UJt dt for $t(u) > 0: 

N E 1 _ 1 
= + — — ) (4.8) 
k— 1 

We note that, almost surely, there exists c > such that 

I^eMI < cr-r, ME, G .4 := {re ie : < r < 00, 101 < \tt} (4.9) 

1 1 M ' 4 



SPECTRAL CHARACTERISATION OF AGEING 



17 



This follows easily from the bound below where u) = a + ib and ./V is any positive integer: 

V 1 f— ^ ifa>0 
1 iw+a!i| — 



IE 



> 



— u + Xj 1 [ {a+x h i)2+b2 ifa<0. 

Let us now introduce the path 7 consisting of the parabolic arcs {— t db it 2 : t > 1} and 
the circular arc of radius v2 around the origin connecting (in anti-clockwise way) — 1 — i 
to — 1 + i. The orientation of 7 is such that — 1 + i comes before + Then, by means of 
(|4.9|) . the inverse formula of Laplace transform and the dominated convergence theorem, 
we get 

f 00 1 

lim F E (x E (t) = Xi) = I e tuJ F(u})du, F{u) := (uj(uj + Xj ) V 

E{-oo h V k=l 10 + Xk 

Note that F(u) is the limit of Fe{cu) as E J, 00, in particular it satisfies (|4.9|) . Moreover, 
F{uj) is the Laplace transform of hmsj,-oo PsCsjsft) = Xj) and trivially 

wF(w) - =^ < c|w|, V|w| < 1 : ueA 

At this point ()4.4I) follows from Proposition IA.1I Moreover, from ()4.4I) and the identity 

N E 

3=1 

it is simple to infer ()4,5I) and (|4.6|) . □ 
4.2. r = e B . 

Note that by choosing to = e^, the random variables X%,X2, ■ ■ ■ introduced in (|4.2|l 
are i.i.d. with distribution p(X)ciX = aX a ~ 1 dX on [0, 1]. Therefore, due to (|4.2|) . we can 
think of Ye(*) as the grand canonical version of the Bouchaud's REM-like trap model. In 
particular, it exhibits the same asymptotic spectral density and the same ageing behaviour: 

Proposition 4.3. For almost all E_, 

1 N E 

lim -77-/ S (e\ =ax a ~ 1 dx weakly in A4([0, 1]). 

EJ.-00 Ne x j 

Proof. By approximating continuous functions on [0, 1] with step functions having rational 
values and jumps at rational points, it is enough to prove that, given < a < b < 1, 

lim -L\{j;l<j< N E , Af> G [a, b]}\ = b a - a a a.s. 

We set 

A E :=\{j- 1 < j < N E , Xj G [o,6]}| = |{i : j > l,e"^ G [e~ E a, e~ E b]}\. 
Then, due to Proposition 12.11 we only need to prove that 

lim = b a — a a a.s. 

El- 00 N E 

To this aim observe that Ne and A E are Poisson variables with expectation (or equivalently 
variance) respectively equals to e~ aE and e~ aE (b a — a a ). Given a positive integer n we 
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set E(n) = ——Inn, i.e. e aE ( n ) = n 2 . It is simple to derive from Chebyshev inequality 
and Boreh-Cantelli lemma that 

1- N E(n) .. 4e(»i) , a a 

lim tttt = 1, nm = o — a a.s. 

raToo e - aE ( n > nToo e ~ aE W 

By a simple argument based on monotonicity, one can extend the first limit to 
lim.Ei-oo -aE = 1 a - s - I n order to extend the second limit to general E we observe 
that, whenever E(n + 1) < E < E(n), 

\A E -A E{n) \ < \{j : e~ E i G [ae^W, ae^^ 1 '] U [be~ E <- n \ be~ E ( n +V] } \. 

Since the r.h.s. is a Poisson variable with expectation of order 0(n), by means of Cheby- 
shev inequality and Borel-Cantelli lemma we obtain 

.. \A E -A E{n) \ 

nm sup , , — = U a.s. 

n ^°° E(n+l)<E<E(n) e {U) 

thus allowing to prove that lim^^oo = b a — a a a.s. □ 
Proposition 4.4. For almost all E_, 

lim n E (t,t w ) = — / - J ° A T — ^ dA ' Vt '^ 

In particular, for almost all E, given 9 > 

lim lim U E (9t w ,t w ) = sm ^ f u - a (l-u) a - l du. (4.11) 

t w ]ooEi-oo 7T J _§_ 

l+s 

Proof. Our starting point is (j4.2j) and the following inequality, valid for any bounded 
function / with E(f(Xi)) = 0: 

kS 2 

P(|AvJ =1 /(X,-)| >S) <2exp(- — — ), W>0,fc = l,2, 



In particular, by conditioning on (see (j4.2j) we get 

P( \Avfzj( Xj )\ >S)< 2exp|-e- aB (l - e" A^)} (4.12) 

It is simple to derive (|4.1U|) from the above estimate, Borel-Cantelli lemma and the integral 
representation 

i r -t A Av Ne e ~ Xjt 



= 5s7 T — i^X dX (4 - 13) 

where 7 is a positive oriented closed path around [0, 1] (see 12. 12^) . Note that the r.h.s. of 
Q4.1UJI corresponds to the function H(t,t w ) introduced in Proposition 12.51 Therefore, the 
conclusion of the proof follows from Propositions 12.51 and 12.81 □ 
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4.3. r | after E [ -oo. 

In this scaling regime, we show that the vague limit of the suitably rescaled spectral 
density is given by the measure ax a ~ 1 dx on [0, oo) and we recover the ageing property of 
the correlation function as before. What is more, however, is that due to the fact that 
we are effectively already at 'infinite times' on the microscopic scale, we get a pure ageing 
function even before taking t and t w to infinity: 

Proposition 4.5. For almost all E, 

N E 

lim lim TqS^S ih =ax a ~ 1 dx vaguely in -M([0, oo)). 
to 10 El- oo f \ 

Proposition 4.6. For almost all energy landscape E, given positive t,t w , 

lim U E (t,t w ) = — — J ~ X \ X dX (4.14) 

and 

lim lim U E (t,t w ) = sin ( 7ra ) „-«(! -u)"" 1 ^ ^ere = — . (4.15) 

roi0i?|-oo 7T tu, 

Remark The integral in (|4.14|) exists due to Lemma 14.81 

Due to Proposition 12.11 Proposition 14.51 follows if one is able to prove that Tq Ylj=i 
converges vaguely to ax a ~ 1 dx on [0, oo) when taking the (ordered) limits E [ — oo, tq { 0. 
This is the content of Lemma 14.71 below concerning the self-average property of Poisson 
point processes with finite intensity measure (compare it with Lemma 4.16 in |3j). Finally, 
the proof of Proposition 14.61 is based on (and given after) the technical Lemmat a 14 . 71 POl 

Lemma 4.7. Let M > and let f be a bounded continuous function on [0, M]. Then 
there exists S > such that for almost all energy landscape E 

VS Yl f( x *) ~ / f(x)<*x a - 1 dx\ < cr 5 , Vr > (4.16) 

Xi<M ^° 

for a suitable positive constant c. 

Proof. Let Xi, X2, ■ ■ ■ and i%m be as in (|4.1|1 . Due to (j4.1|) and since Var(nM) = E(riAf) = 
(M/T Q ) a , 

\{j : Xj < M }\ 
{M/t q Y 

In particular, given 7, s > with 2s — 7a < —1, by means of Borel-Cantelli lemma we 
obtain that for almost all energy landscape E there exists c > such that 

10 ■ < M}\ 



{M/T ) a 

Due to the above estimate, 



< ck s , VA; = 1, 2, ... where tq := k 7 . 



\ro Yl ~ M a Av Xi < M f( Xl )\ < cfc-'H/lloo, V/c G 1,2,... where r := k^ 

(4.17) 



20 



ANTON BOVIER AND ALESSANDRA FAGGIONATO 



where Av Xi <M denotes the average over the set {xi < M}. As done for (|4.12j) . if < p < 1, 



\ A.v x .< M f(xi) - M' a / f(x)ax a - 1 dx > p) 



<2exp{-(-^) a (l-e- cp )} <2e- cp T ° 



In particular, by Borel-Cantelli lemma, for almost all E, 

*M 



Av Xi < M f(xi) - M~ a / f(x)ax a ~ 1 dx 



o 



< ck~ 



Vfc = 1, 2, . . . where tq := k 



(4.18) 

if s is chosen small enough. At this point (|4.17|) and (|4.18|) imply the assertion of the 
lemma, if To = A;~ 7 for some k = 1,2,... The general case tq > follows easily from the 
uniform continuity of /. □ 

Lemma 4.8. For almost all energy landscapes E, there are positive constants Tq,ci,C2 
having the following properties. If tq < Tq , N > \{j : Xj < 1}| and A 6 (or A = a + ib 
with \b\ < 1 and a > Xj + 1 for all j < N) then 

N 1 

I t oE^aI^ c ii a i~ 2 - ( 4 - 19 ) 

3=1 3 

Moreover, if tq < Tq and M > 1, i/ien 
00 ^ 

r o" E r-TTl ^ C2 l A l a_1 ln ( X + l A ')' 2/ A G 7oo (4.20) 
j =1 \ x j A \ 

X ^ 

T o a V l TT < caM"" 1 InM, i/ A G 7oo or K(A) = M + 1 (4.21) 

^— ' x 7 - — A 

ay<M 1 J 1 

Proof. It is convenient to introduce the non rescaled Ppp ^2 i 5 Vi , with yi := e~ Ei , having 
intensity measure ay a ~ l dy on (0, 00). Moreover, we fix here (3 > 2 and < 7 < (3/2 — 1 
and we define 

N n : = j {j : vP/ a < Vj <{n+ l)M a } \. 

for n positive integer. Then a simple application of Borel-Cantelli lemma implies that, 
for almost all energy landscapes E, 

7—^ --I <cn~\ Vn = l,2,... (4.22) 

(n + 1)p — nP 

In fact, by Chebyshev inequality and since N n is a Poisson variable with Var{N n ) = 
E(N n ) = {n + lf-nP, 

F(\N n /E(N n ) - 1| > n~ 7 ) < cn 27+1 ^ 

Moreover, by Borel-Cantelli lemma and a simple argument based on monotonicity, it 
is simple to prove that there exists 5 > such that, for almost all E, 



\{j ■ Vj < u}\ _ 



< KU~ 5 Vll > 1. 



(4.23) 



for a suitable positive constant k. 
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In what follows we write A = a + ib. Then 

N N 



1 



\T- 

j=i j 



E 



+ 



^ (xi - a) 2 + b 2 J \^ (xj - a) 2 + b 2 



N 

(E 



(4.24) 



In order to prove (|4.19|) we assume (|4.22|) and ()4.23j) to be valid and let < r < Tq < 1 
where Tq is such that u a — ku~ 7 > for all u > 1/tq. In particular, {xj : Xj < 1} ^ 0. 
By dOH) , if N > \{j : Xj < 1}| and A G 7oo 



70 



N 

IE 

i=i 



X, — A 



> 



'tf £,,<! (xj J%* >c|A|- 2 r "|{x 3 - : x, < 1}| if |6| > | 

if H < i 



r i? Sxj<i - ct q\{ x j '■ x i - !}l 



At this point, (|4.19j) for A S 7oo follows from (|4.23j) . The case \ = a + ib, with |6| < 1 
and a > xj + 1 for all j < N, can be treated similarly. 

It is simple to derive (j4.2U[) and (|4.21|) from estimates (|4.25|) . . . . , ()4.31[) below valid for 
almost all E: 



if a < 100, 1 < M and A = o + ib G 7oo then 



^Ej^^. 



o-l. 



if tq is small enough and a > 100 then 



< ca 



T » E ^ 

r Q V 1 — - — -<ca a ~ x \ua 

u U . _ n \ — 



T, 



\Xj — a\ 

I, ^ i — rr < ca" -1 if A = a + ib G 7 C 

— A 

a-l<A<a+l 1 J 1 

o V 1 " r< carina, 

u ^ \xi - a \~ 

a+l<Xj<2a ' J ' 

- Q V -. — - — r < ca-^M"" 1 , if M > 2a 



Xj>M 1 J 



Let A G 7oo with a < 100. Then, due to flQ5j) , 

t »° E^-tai £ "»° I{j ' : w-^: )l - c ' 



7"0 



while, due to (jQ3Jl . 



^^^<E^=r E 



n^ 1 ^ < c" 



(4.25) 
(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 



(4.32) 



22 



ANTON BOVIER AND ALESSANDRA FAGGIONATO 



thus proving (|4.25j) . The proof of (|4.26j) follows the same arguments of (j4.32j) . (|4.27|) is 
a simple consequence of Q4.23JI . The l.h.s. of (|4.29|) can be bounded by T^\{j : a — 1 < 
ToUj < a + 1}| and (|4~2"2"j) allows to conclude the proof of (|4~2"$|) . 

The proof of (j4.28j) . (|4.3()j) and (|4.31|) can be easily derived from the following estimate. 
Let 1 < A < B with B<a-lorA>a + l, then (j4~22j) implies 



\ '* ^ In I x 



A<Xi<B 1 1 n=n._ 1 u 1 ,/u 1 u 1 

= c «" / - — — — 7r, — -dy 

where n_ = LC^Ao)"^ 3 ] — 1) n + = L(-^/ r o) a ^J + 1, u = ^- — 1, f = n + + 1 (we assume 
To small enough in order to exclude the singular point in the above intervals of sum and 
integral) . □ 

Proof of Proposition UTTH In order to avoid confusion we underline here the dependence 
on to = e E ° by writing HE,E (t,t w ) in place of IL^i, t w ). Our starting point is given by 
the integral representation ()2.12j) : 



1 f e twA r o 2^=1 37 



-A 



Il E , Eo (t,t w ) = — / — ' \ dX. (4.33) 



Let us choose E_ satisfying Lemma [4 .81 Then, due to the exponential decaying factor e~ twX 
and to Lemma l4.8l if ro < Tq and E is small enough such that T$e~ E > 1, the path integral 
7^ in (|4.33|) can be substituted with 700. At this point, ()4.14j) follows from Lemma |4"%1 
and the Dominated Convergence Theorem. 

In order to prove ()4.15j) . given a positive integer M, we set 

^ ^ — Xjt 

If e~ twX t oz2 Xj <m ^px 
9M,E {t,t w ) := — / — — — dX 

J L M /-jXj<M Xj—X 

where Tm is the positive oriented path having support 

supp(r M ) = {A G C : |A - x\ = 1 for some x G [0, M]}. 

Then, by applying Lemma 14.81 whenever ro < Tq 

I lim Ue e (t, t w ) — gM,E (t, t w )\ < cM " 1 In M, VM G N+. 
Ei-00 

Let us assume that E_ satisfies (|4.16p for all M G N+ and for f(x) = ^zr\ 01 fi x ) = ' x ~^\> 
for all A in a countable dense set of Tm and for all rational positive t. Then, by a chaining 
argument, we get 

lim gM,E (t,t w ) = gM{t,t w ), Vt,t w > 

Eol-00 

1 f e-*» A C ^x a ~ l dx 

9M(t,t w ) = — - J M A - — dX, (4.34) 

2m JTm A Jo j^x^dx 



where 



thus implying 



limsupl lim U E>Eo {t,t w ) - g M (t,t w )\ < M a_1 lnM, VM G N 4 
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At this point let us observe (see the proof of Lemma l4.8j) that there exist c, d > such 
that for all M £ N+: 



rM 

\ 


C¥ — 

X 


i 


Jo 


A - 


— > 

X 


rM 


x a - 


-1 


Jo 


|A- 


X 


POO 


x a ~ 


-1 


Jm 


|A- 


x\ 



ds\>c\X\ , VA G Tm U 7oo 

T cte<c', VAGr A/ U7oo (4.35) 



da; < c M a ~ InM, VA E 7oo . 

From the above estimates we infer 

\g M (t,t w ) - g(t,t w )\ < cM^bM (4.36) 

where 

Using the analytic properties of the integrand in the r.h.s. of (|4.37j) . one can show that 
g(t,t w ) = g(t/t w , 1). In order to compute g(6, 1), we observe that for a suitable positive 
constant c \g(9s,s) — ghi{6s,s)\ < cM a-1 mM, for any s > 1 (in fact, the constant c in 
()4.36|) can be chosen uniformly if t w > 1). By the results of Subsection 12 . II (compare ()4.34j) 
with LT(t, t w ) in Proposition 12. 5() we get 



lim ghiids, s) = r.h.s. of ()4.1 

si oo 

thus concluding the proof. □ 

5. Other correlation function when t j after E [ — oo 

As stated in Proposition I4.6| the standard time-time correlation function H E (t,t w ) has 
trivial behaviour after taking the limits E { —oo, tq j 0. For physical reasons, it is more 
natural to disregard jumps into states with physical waiting time T\ = e Ei much smaller 
than ro, since we assume that the physical instruments in the laboratory have sensibility 
up to the time unit ro. Therefore, let us fix 5 > and consider here the more natural 
time-time correlation function 

U^(t,t w )=F E (x E (u)>5 VuG (t w ,t w + t\: x E (u) ^ x E {u~) ) . 
The main result of this section if the following one: 
Proposition 5.1. For almost all E_ 

lim sup lim lim I Hit (t, t w ) — H E (t, t w ) I = 0. (5.1) 

T°o f>o S0J.-00 El— 00 

In particular, for almost all E, 

lim lim lim Ilffi (0t w , t w ) = sm ^ a "> f u - a (l - u^alu, V# > 0. (5.2) 

t^TooEol-ooEI-oo V ' 7T J _§_ 

1+9 

Note that the correlation function Tl^\t,t w ) is the analogous of Yl^j\t,t w ) of Propo- 
sition As for the proof of Proposition 13.11 a useful observation is that, given 5 > 0, 

limlim lim F E (x E (t) > S) = a.s.. (5.3) 

tlooT lOEl-oo V ' 
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We can prove a stronger result concerning to the phenomenon that with high probability 
the system visits deeper and deeper traps. In fact, note that by (|2.13j) 



p, (x , w >o=^/ Vf£^ 

lb L~tj=\ Xj—\ 



Then, by reasoning as in the proof of Proposition 14.81 and using the results of Section \2. 21 
one can easily show the analogous of Proposition 12.91 

Proposition 5.2. For almost all energy landscapes E_, 

limlim lim t l ~ a ¥ E { x E (t) > 5) = (5.4) 
iToo to|0 B|-oo v ' c(a) 

where 

f°° x a ~~^dx r°° 
B ( 5 ) ■= foo^-i , . c («) : = / iT^dy 
Jo — dx J o 

Proof of Proposition IBTT1 Trivially, (|5.2j) follows from (|5.1j) and Proposition I4.fil Since 
Eq I — oo after E { — oo, we assume that E < Eq and define 

D E:Eo :={i :E<Ei< E }. 

This set corresponds to the small traps where we allow the particle to jump in. By l|5.M|) 
and the same arguments used in the proof of Proposition 13,11 for i = 1 (with exclusion of 

the last step since here lirngj-oo ^ff^ = 1 a.s.) it is simple to derive the assertion of the 
proposition from Lemma 15.31 □ 

Lemma 5.3. For almost all energy landscapes E there exist positive constants p,c (in- 
dependent of E,Eq) satisfying the following property. Whenever \{i : Ei > Eq}\ > 0, 

limsup-^- V (p E ,Ea(i,t) < ce~ pt , Vt > 0, (5.5) 

where D EiEo '■= {i '■ E < Ei < Eq}, for E < Eq, and the function <Pe,e is defined as 

<pE,E (i,u) := F E (Y E (u) G D E , Eo G [0,s] \ Y E (0) = i), (5.6) 

Proof. Let us assume that E < Eq, \{i : Ei > Eq}\ > and, without loss of generality, 
<5 = 1. 

We fix t > such that e~ ai < I and define 

W 1)E :={i : E<Ei<E + £}, N 1)E := \W 1)E \ 
W 2>E :={i : E + i<Ei< Eq}, N 2jE := \W 2 , E \ 

Note that D E<Eo = W\ tE U W 2 , E and that N E , N± iE , N 2 , E are Poisson variables having 
expectations e~ aE , e~ aE (l — e~ ai ) and e~ aE ~ ai — e~ aE ° . In particular, for almost all E_, 

r N E 1 N 1>E _^ N 2tE _ ae 1 

hm ? = 1, pi := lim — — = 1 — e , p 2 := hm = e < - 

Ei-oo e~ aE ' F Ei-oo N E F Ei-oo N E 2 

(5.7) 

We observe that h := N E — N\ >E — N 2>E is a positive integer independent of E and 
Xi > e E °~ E ~ e if i € W\ E , while Xi > 1 if i € W\ E . Let us introduce a new random walk 
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Y E {t) on S E whose infinitesimal generator 1L* E is denned as with Xj replaced by x* 
defined as 

Xi tfigD E ,E , 
A E :=e E °~ E ~ e i£i€Wi tE , 
1 if i G W 2 ,e 

We denote by P* E the probability on path space associated to Y E (t) when having initial 
uniform distribution and we set 

<p EjEo (i,u) := P E (Y E (u) G D EtEo G [0,s] | Fj(0) = i). (5.8) 

By a simple coupling argument, one gets ip EjEo (i,u) < ip* E E (i,u). In particular 

^ £ ^,Bb(*.*)<*:= ^ E ¥>WM)> V»eD B) Eb.V*>0. (5.9) 

At this point it remains to estimate <3?. In order to simplify notation we write simply 
D, N, N\,N 2 ,A, by dropping the index E. Moreover, we consider the following realization 
of the dynamics of Y E : after arriving at a site i, the system waits an exponential time of 
parameter x* and the it jumps to a point of S E with uniform probability. In particular, 
jumps can be degenerate, i.e. initial and final sites can coincide. 

We claim that 

$ = + $ 2 + $ 3 

where 



(5.10) 

fc 1= 0fc 2 =l V 7 U ' " 



* 3 := §«p{-*(l - *)} + $«p{-t(l - f )} (5.12) 
The above identities can be derived from the probabilistic interpretation of k±, k 2 as 

ki = |{ jumps performed before time t having starting point in Wi}\ 
and from the following simple identities: 

n— 1 

P(Ti + T 2 H + T n G [z, * + dz) ) = e- K ^ n - — dz 

v ' (n — ly. 



7 

k z ,.n 



T\ + T 2 H 1- T n < z and T\ + T 2 H hT n + T n+ i > z) = e " " /v 



rt! 

where z > and T\,T 2 , . . . are independent exponential variables with parameter k. 
Finally, we only need to prove that, for suitable positive constant c,p > 0, 

limsup$i < ce~ pt , Vi = 1,2,3. (5.13) 

E|-oo 

We give the proof in the case i = 1, the case i = 2 is completely similar while the case 
i = 3 follows directly from (|5.7|) , 
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We fix 7 : a < 7 < 1, set fco := A 1 and write 

where is the contribution to <I>i of addenda in the r.h.s. of (|5.1U|) with 1 < k% < kq 

and k 2 > 0. 

If k% > fco then 

N 1 \ki (N-N 2 \ k ^(_, n \ki < ,N - N 2 \ kl (, 1 \ A ~' 



NJ V N ) V N-N 2 ) ~\ N ) \ iV 
thus implying that 

*>*<^^,tf 2 )(l-i) <(l-_) 10 as^i-oo. 

where $i f ^fi , JV 2 ) is defined as in the r.h.s. of with N replaced by N - N±. 

Note that is does not exceed 1 since it corresponds to the probability of a certain event. 
Let us now consider the term &f k °. To this aim, since ( kl k ^ 2 ) < 2 kl+k2 , 

^<k < V ^ ( 2AN, k 2tN 2 k2 I(Q,t) 
1 - ^ ^ 1 iV j 1 N ' (h - l)\k 2 \ 

where 

I( Wl ,w 2 ):= e-^- u) - Au u kx - l du. 



Fix m > with 2p 2 + 2mpi < 1 (recall that 2p 2 < 1). Then trivially 

/(0,^) (^)" (5.14) 

From such a bound, one gets immediately 

EEC— ) (~ ) _ 1)!fc2 , ^(0, -X) < c exp|-t(l - - - 2-m - 2-j ) 

(5.15) 

The last expression, when E [ — oo, converges to cexp(— t(l — 2p\m — 2p 2 )) , in agreement 
with fl5~T3]) . 

In order to estimate the integral I(^,t) = e~ l fLt e~( A ~ 1 ) u u kl ~ 1 , we observe that 

A p-zs _ p -zw p -zs i p —zw i 

e- zu u n du = (-l) n -r- ~ < (s + -T + ~ (w + -)", Vs, w, z > 0, 

dz n z z z z z 

thus implying the bound 

I^T,t) < ce-'A-^imt + l)* 1 " 1 + e" A *(A - l)" 1 ^ + -y^—)^ 1 

.ri ./I J. 

The contribution of ce - '^ 1 (mt + l)^ 1 " 1 to $f fc ° can be treated by means of estimates 
similar to the ones leading to (|5.15|) . 

In order to conclude we only need to show that 

-At f 2AtNi,ki-l/2tN 2 ,k 2 1 | n t-i I 

e EEHH (— ) - l) !fc2 . + °» »^i-oo. (5.16) 

fcl=l /C2 = l 



SPECTRAL CHARACTERISATION OF AGEING 



27 



To this aim observe that 

AT 



r.h.s. of GDI < e- At+2t % £ (^T 1 )" 1 



fei=0 

< c^e^A^^tA)^ = c(t) exp {- At + j In A + A 7 ln(2L4)} 
Since < 7 < 1, we get (|5.16[) . thus concluding the proof. □ 

Appendix A. Laplace Transform 

Proposition A.l. Let G(t) be a bounded measurable function on (0,oo) and let us con- 
sider the Laplace transform 

POD 

G(uj) = / G(t)e- tuj dt 
Jo 

well defined if ?R.(lj) > 0. Let us define 

A := {re ie : < r < 00, |0| < ^vr} (A.l) 

Suppose that G can be analytically continued to C \ (—00, 0] and that there are positive 
constants 7, (3, a, c and B6l such that 

\G{lj)\ < c|^r 7 Muj € A, \u\ > 1, (A.2) 

\uPG(w) - B\ < c\u\ a Vlo 6 A, \u>\ < 1. (A.3) 

Then, 

lim s l -^G(s) = — — where c(8) := / y^e^dy. 
sToo c(/3) Jq 

Proof. If we set H(s) := ^py s/3_1 with s > 0, then the Laplace transform H(u>) is well 

defined for 3^(cj) > 0, H(uj) = Blu^P and trivially H(uj) can be analytically continued to 
C\(-oo,0]. 

By the inverse formula for Laplace transform (see Chapter 4, Section 4 in |16j). we have 

-1 px+iK 

G{s) = lim / e SUJ G(uj)duj, Vs > 0, x > 0, (A.4) 

K^oo 2m J x _ iK 

where oj runs over the vertical path connecting x — iK and x + iK. The above formula 
remains true if substituting G with H. Therefore, 

s 1 " /3 G(s) - — — = lim - / e s "(G{uj) - H{uj))oIuj, Vs > 0, x > 0, (A.5) 

c(P) K^oo 2m J x _ iK 

Let p := min(7, (3)/2. Fix a positive number x and, given K and s, define the following 
paths (see figure below). 

is the vertical path from x — iK to x + i-fT. 71 + is the segment from — s _1 + s~ l i to 
— 1+1 72,+ is an arc from —1+i to —K p +iK given by the parametrisation z{t) = —t-\-it l l p 
with t S [1, K p \. 73 j+ is the horizontal segment from —K p + iK to x + iK. For i = 1, 2, 3 we 
define the path 7^ _ by considering the reflection of 7^4. w.r.t. the real axis and inverting 
the orientation. Let 70 be the positive-oriented circular arc of radius s _1 from — — s~ 1 i 
to —s" 1 + s~ 1 i crossing the axis of positive real numbers. 
Note that the above paths depend on s and/or K. 
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-KP + iK 




x + iK 



IK 



-KP - iK 



x — iK 



Figure 1. The integration paths 



Because of analyticity, the integral over 7^ of e su G(w) is equal to the sum of the inte- 
grals over 73,-, 72,-, 7i,-> 70j 7l,+> 72,+, 73,+- The same is valid with G replaced with H. 



By (|A.2j) we have that 



\e su G(uj)\ \dw\ < ce sx K p -^ [ as K | 00 



(A.6) 



73, ± 



and, for a suitable rational function /, 



,1-/9 



s * / \e SUJ G(tu)\ \du\ <a* " I \r 

'72,± J I 

f 00 
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1-/3 / L^-t+it 1 /") 



G(-t + + p~ x rf%) \dt 



/CO 
e- st f(t)dt < c's 1 -^-^ I 0, as s T 00. 



(A.7) 



Similarly, it can be proved that the corresponding integrals with G substituted with H go 

to by taking the limits K f 00, s j 00. 

Let us now estimate s 1- ^ J _i e~ st t a ~ l3 dt by dividing the path of integration in two 

paths. Choosing < 5 < 1, 

. a -i+i 
g i-/3 / 



r „t t a~f3 dt < cs l-f3\ s (-l+5)(l+a-P) _ g -(l+a-/3)\ < J g -a + J s -a+S(l+a-f3) 



,1-/3 



-l+«5 
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for a suitable rational function g(s). In particular, choosing 5 small enough, the above 
upper bounds imply 

lims 1 ^ f e- st t a ~ f) dt = 0. 

This result, together with ()A,3|) . implies 

1 



lims 1 "^ / e SUJ (G(Lo) -H(u))dw = 0. 



71, ± 



Trivially, by ifOjl . 



4-/3 1 / p sw t 
'To 



> -| / e sa; (G(u;) - < s" a | as s | oo. 



The Proposition now follows from the estimates above and ()A.5|) , 



□ 



Appendix B. Perturbation Theory 

In this appendix we comment on a paper by Melin and Butaud L 22J where the eigenvalues 
and eigenfunctions of the generator of our model were computed using perturbation theory. 
As pointed out earlier, these results are at variance with our exact results, and it may be 
worthwhile to point out the flaw in their arguments. Melin and Butaud write the generator 
L defined in as L = T + j^T^ where 



T := 



(xi 
x 2 







o \ 



x N J 



( 



r (i) := 



-X 2 



\-x N 



-Xl 

-x N 



-x\ \ 

-xi 

-XN J 



(B.l) 



The factor 1 /N in front of the second term encourages them to consider this term as a small 
perturbation. Both T and are symmetric operators on L 2 (/j,) where := xj 1 . We 
denote < •, • > the scalar product in L 2 (n) and assume x±, . . . ,xn to be distinct positive 
numbers. 

Given an operator A : L 2 ([i) — > L 2 (,//) we write \\A\\ for its operator norm. Because of 
symmetry, ||T|| and ||TM[| are given by the maximum of |A|, with A eigenvalue. Trivially, 
T has eigenvalues x±, . . . , xjy and T(ej) = xiei where ei, . . . , ejv is the canonical basis of 
M. N , while has eigenvalues 0, — {x\ + x 2 + • • • + xn)- 

Given z G C we can define the holomorphic function T(z) = T + zT^\ A natural 
condition in order to apply perturbation theory to T(z) (see [231 > chapter II) is 



\z\ < 



d_ 

2a 



(B.2) 



where 



d = inf \xi 



X 



3\> 



a :- 



min||T (1) 



xi + x 2 H h xn 
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In this case, we can conclude that T(z) = T + zT^ has N eigenvalues \i(z), . . . , Xn(z) 
with A fc (z) = E"=o4 n) ^. l4 n) l < ^{2/dY- 1 and 

m _ < T^e k ,e k > 



K: 



< e k ,e k > 



< e fc ,e fc >< e 3 -,ej > 



Similar series exist for the perturbed eigenvectors. 

However, the crucial condition (|B,2|) is hardly satisfied when z = jt, since it reads 

Av^ =l Xj < inf \xi — Xj\ (B.3) 

while a.s. the l.h.s. of (|B.3|) has non zero limit and the r.h.s. converges to like l/N. 

The fact that the conditions for the application of perturbation theory are violated 
explains why its predictions are incorrect. This discrepancy happens not to be too obvious 
as far as the eigenvalues are concerned (which are caught between the diagonal elements of 
the generator and thus are somewhat similar to them, but the shape of the eigenfunctions 
is sharply different). 

Namely, by of Proposition 12.11 when j ^ 1 and Xj-i,Xj are very near each other, the 
eigenvector ipv) related to the eigenvalue A,- : Xj-± < Xj < xj has two main peaks of 

opposite sign given by ijjj-i an d , this is very different from the predictions of [2^] (see 
their Figure 4). 

Appendix C. Complex integral representation 

Let L be a Markov generator on the state space S := {1, 2, . . . , N}, reversible w.r.t. a 
positive measure /i. We can think of L as a linear operator on W N , symmetric w.r.t. to 
the scalar product (•, -) M where 

N 

(o,6) M = y^p,{i)aibi 
i=i 

In what follows we endow with the scalar product (•, ■)„ (and not with the standard 
Euclidean scalar product). Since L is symmetric, we can orthogonally decompose M> N as 
R N = Wi W 2 © ■ ■ ■ © W m such that L = YlT=i ^kP\v k , where P Wk denotes the orthogonal 
projection on W k and Aj / Xj if i ^ j. Given A G C \ {Ai, . . . , A m }, we write -R(A) for the 
resolvent 



k=l 

Then, the Residue Theorem implies the integral representation 



e~ tL = / e- tx R(X)dX (C.l) 



2iri 



7 



where 7 is a positive oriented loop containing in its interior Ai, A2, • • • , X m . 

Given a probability measure v on S, we denote by ¥ u the probability measure on the 
path space associated to the continuous-time random walk Y(t) on S with generator L 
and initial distribution v. Fix j G S and let v G M N be such that V{ = 5ij. We write 
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for the Radon derivate, i.e. = ^y- Then the symmetry of L w.r.t. the scalar 

product (v)^ implies 

p,(y(t)=i) = E^)( e - iL ) fe!i = M| !e -N =^%v) 

u(k) 



k=l 

N 



By plugging (|C1|) in the r.h.s. of (|C2|) we get the integral representation 

N 

>(Y(t)=j) 



N N OA r 

7 — 1 k — 1 



(C.2) 



In particular, given h function on S 

N N 

I ;/[/,■!/ 

,(A)dA 

thus allowing to get an integral representation for H(t, t w ) := Pj,( no jump in [t w , t w + 
If we set ^(i) = Tj = x" 1 and = N -1 (uniform initial probability), then 

F v (Y(t)=j) = ^J ^{^E^^W}^- ( C - 4 ) 



Let us consider now the special case given by the Bouchaud's REM-like trap model 
where L := L^r is defined in ()2.3j) and v is the uniform distribution on S. Note the all the 
integral formulas obtained in Section [21 can be derived from the following one: 

P - (r " )=J '-2k/ T e " A (A-^(A) rfA < C5 > 

where <p(X) = Y2k=i \~x • ^ n w hat follows we prove that ()C.5|) corresponds to (|C.4|) , 

We know already that det(AI — L) has distinct zeros given by the N distinct zeros of 
0(A). In particular, it must be 

det(AI — L) = I^(A) l[(X - x 3 ) = 1a £ [J (A - x,) (C.6) 

Given a matrix ^4 we write L4]i,j for the determinant of the matrix obtain from A by 
erasing the i-th row and the j- column. Since 

R ^ k{X) ~ det(AI-L) 
and due to (|C.6|) . in order to derive (|C.5|) from (|C4|) we only have to show that 

£(-l)i+*+^[AI-L] w = J] (A -a.) (C.7) 

In order to prove the above identity observe that [AI — L]/% j is a polynomial of degree N — 1 
if k = j, otherwise it has degree N — 2. The l.h.s. of (|C.7|) is a monomic polynomial of 
degree N — 1. At this point we only have to verify that x s , s ^ j, are zeros of the l.h.s. 
of (|C.7j) . This is trivial if one observes that the l.h.s. of (|G7j) is the determinant of the 
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matrix obtain from AI — L by replacing the j-column with the vector w with Wi = — for 

Xj 

i = 1, 2, . . . , N. It is simple to verify that, if A = x s for some s j, the j—th row and the 
s—th row in such a matrix are proportional, thus implying the thesis. 

Acknowledgements. The authors thank J. Cerny and V. Gayrard for useful discussions. 
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